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Abstract 

We study the notion of structured realizability for linear systems defined over graphs. 
A stabilizable and detectable realization is structured if the state-space matrices inherit 
the sparsity pattern of the adjacency matrix of the associated graph. In this paper, we 
demonstrate that not every structured transfer matrix has a structured realization and 
we reveal the practical meaning of this fact. We also uncover a close connection between 
the structured realizability of a plant and whether the plant can be stabilized by a 
structured controller. In particular, we show that a structured stabilizing controller 
can only exist when the plant admits a structured realization. Finally, we give a 
parameterization of all structured stabilizing controllers and show that they always 
have structured realizations. 



1 Introduction 



Linear time-invariant systems are typically represented using transfer functions or state- 
space realizations. The relationship between these representations is well understood; 
every transfer function has a minimal state-space realization and one can easily move 
between representations depending on the need. 

In this paper, we address the question of whether this relationship between representa- 
tions still holds for systems defined over directed graphs. Consider the simple two-node 
graph of Figure 1. For i — 1,2, the node i represents a subsystem with inputs Ui and 
outputs yi, and the edge means that subsystem 1 can influence subsystem 2 but not vice 
versa. The transfer matrix for any such a system has the sparsity pattern <Si. 

For a state-space realization to make sense in this context, every state should be as- 
sociated with a subsystem, which means that the state should be computable using the 
information available to that subsystem. In the case of S\, the states that determine y\ 
must only depend on u%, while the states that determine 2/2 may depend on both u\ and 
U2- Consider the following example, which belongs to the graph constraint of Figure 1. 
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Figure 1: A two- node graph and its adjacency matrix. 



More generally, we say a system is S -realizable if it has a stabilizable and detectable 
realization for which A, B, C, D each have the sparsity of S. A more formal definition is 
given in Section 3.3. 
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While this definition seems natural and one might expect every structured transfer 
function to have a corresponding 5-realization, it is not the case in general. For example, 
Gi defined in (2) belongs to the graph constraint of Figure 2, but no ^-realization exists. 
A proof is given in Appendix A. This example shows that there is no immediate relation 
between the sparsity of a transfer matrix and that of its state-space matrices. 
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Figure 2: A four-node graph and its adjacency matrix. 



Our main result is that <S-realizability is necessary for 5-stabilizability; finding a stabi- 
lizing controller that has the same 5-structure as the plant. Non-realizable systems exist, 
but such systems cannot stabilize or be stabilized by other structured systems. 

The paper is organized as follows. In Section 2, we cover some related work in the 
literature and in Section 3 we cover basic definitions and concepts touching on systems 
over graphs, realizability, and stabilizability. Our main results are in Sections 4 and 5, 
which are followed by concluding remarks in Section 6. 



2 Literature Review 

There is a large body of work exploring control systems defined over graphs. For a broad 
class of systems, synthesizing an optimal controller can be reduced to solving a convex 
optimization problem [5, 7, 9]. 

Several solution approaches have been reported. LMI methods [1, 8, 12, 20] work 
directly with state-space realizations for the plant and controller. Vectorization [10] avoids 
the sparsity constraint by reshaping the transfer matrix. Alternatively, one can solve a 
sequence of finite-dimensional convex problems whose solutions converge to the optimal 
structured controllers [7]. 

Surprisingly, the issue of structured realizability is not addressed in any of the afore- 
mentioned works. Indeed, the LMI methods assume that structured realizations for the 
plant and controller always exist, while transfer function methods make no mention of 
state-space — so it's conceivable that a transfer function method might generate a con- 
troller which has no structured realization! 

In this paper, we consider directed graphs with delay-free links. Our framework is simi- 
lar to that of Swigart and Lall [15]. For many such problems, explicit state-space solutions 
have been directly computed [3, 4, 6, 13, 14]. In these works, a structured realization is 
assumed for the plant, and the optimal controller turns out to have a structured realiza- 
tion as well. The topic of whether such realizations should exist in general, or whether 
suboptimal structured stabilizing controllers should always have structured realizations 
is not discussed. 

To the best of our knowledge, the only existing work dealing with structured realizabil- 
ity is the work of Vamsi and Elia [16, 17]. In these papers, the authors provide sufficient 
conditions for structured realizability, as well as an LMI approach for controller synthesis 
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that guarantees realizability. However, they give no example of a provably non-realizable 
system, and the sufficient conditions provided are potentially very restrictive. 

In Section 1, we showed that structured realizability is a meaningful concept by giving 
an example of a non-realizable transfer function. In the sections that follow, we provide 
some very general results; a parameterization of all structured stabilizing controllers, and 
a proof that 5-realizability is necessary for 5-stabilizability. 

3 Preliminaries 
3.1 Directed Graphs 

A directed graph is a set of nodes V = {1, . . . , N} and edges £ C V x V. If (i, j) G £ , we 
say that there is an edge from i to j, and we write i — >• j. We make several assumptions 
regarding the graph. 

Al) Self-loops: for alH G V, i -> i. 

A2) Transitive closure: if i — > j and j — > k, then i — V k. 

A3) There are no directed cycles of length 2 or greater. 

These assumptions are natural, and will be further justified in Section 3.2. Given a graph 
satisfying the assumptions above, define the adjacency matrix S G {0, l} WxlV 



Under assumptions ^41-^43, one can always relabel the nodes such that S is lower- 
triangular with a full diagonal. Therefore, we will assume a lower-triangular S from 
now on. See Figure 2 for an example of a graph and its associated adjacency matrix. 

3.2 Systems over Graphs 

The systems considered in this paper are linear, time-invariant, continuous-time, and 
rational. We denote the set of proper rational transfer functions as 1Z P . If all the poles 
of G G 1Z P are contained in the open left-half plane, G is stable, and we write G G TZT-Loq. 

For systems defined over graphs, additional notation is needed to keep track of the 
input and output partitions. 

Definition 1 (Index sets). Suppose we have a graph S G {0, l} JVxAr . An index set k is 
a tuple (fci, . . . , fcjy) of nonnegative integers. We also define the set of nonempty indices 
as Q k = {i G V | ki ^ 0}. 

Definition 2. Suppose we have a graph S G {0, l} NxN with associated index sets k and 
ra, and F is a commutative ring. We write A G S(F, k,m) to mean that 



If j — > i, then Aij G F kiXm i . Otherwise, A^ =0. If k and m are clear by context, we 
simply write A G S(F). 




otherwise 
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A plant G denned over the graph (V, £) is written in its most general form as G G 
S(1Z p ,k,m). Intuitively, this means that if i — > j, then input affects output yj, and 
the associated transfer function is Gij. If there is no edge from i to j, then dj = 0. 

We seek controllers that obey the same structure. Namely, if i — V j then the control 
signal Uj may depend on the measurement y i: and the associated transfer function is K^. 
If there is no edge from i to j then = 0. 

We can now see why Assumptions ^42-^43 make sense. The graph (V, £) represents 
information flow; if a subcontroller i shares what it knows with subcontroller j along the 
link i — > j and similarly for j — > k, one would expect subcontroller k to have access to the 
information from subcontroller i as well. Directed cycles can also be removed by treating 
all nodes involved as a single node. 

Remark 3. Assumption A2 leads to a useful algebraic property. If G\ G S(R pi k, m) and 
G2 G S(lZ p ,m,p) then G1G2 G S(TZ p ,k,p). In other words, structured transfer functions 
as defined above form an algebra. 

3.3 Structured Realizability 

Structured realizability is a core concept in this paper. Roughly speaking, we are inter- 
ested in finding a state-space realization for a structured plant comprised of matrices that 
are also structured. 

Definition 4 (Structured realizability). A transfer function G G S{lZ Pl k,m) is said to 
be S(k, m) -realizable if there exists an index set n and matrices A 6 S(M.,n,n), B G 
<S(R, n, m), C G 5(K, k, n), and D G 5(M, k, m) such that (A, B,C, D) is a stabilizable 
and detectable realization for G. When m, and k are clear by context, we will simply 
write that G is S-realizable. 

Note that we allow the index set n to have zero-entries. If a component of n is zero, it 
means that no state is associated with that subsystem, and the corresponding rows and 
columns of A collapse. 

A key aspect of this definition is the requirement of stabilizability and detectability. 
Indeed, we will see in Section 4.1 that one can always construct a trivial realization with 
the correct structure by introducing duplicate (non-minimal) states. 

3.4 Stabilizability 

Throughout this paper, we use the conventional notion of internal stability for feedback 
interconnections [2, 21]. If a plant G and controller K are connected in feedback as in 
Figure 3, then K is stabilizing if and only if the map {u\,u\) M> (2/1,2/2) is proper and 
stable. When dealing with decentralized systems, there are other natural ways to define 
stability. For a further discussion on this topic, see Section 6. We now state the formal 
input-output definition of stabilization. 





Figure 3: Feedback loop with inputs and outputs added to the feedback path. 
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Definition 5 (Stabilization). Suppose G G Tl kxm and K G 1l™ xk . We say that K 
stabilizes G if 

(i) I — G(oo)K(oo) is invertible, and 



(n) 



I -G 
-K I 



There is also a useful state-space characterization of stabilization, which we state as a 
proposition. 

Proposition 6. Suppose G G 1Z kxm , K G TZ p nxk have realizations given by (A, B,C, D) 
and (Ak,Bk,Ck,Dk) respectively. Then the following are equivalent. 

(i) (C,A,B) and (Ck,A k ,B k ) are stabilizable and detectable, and K stabilizes G. 

(ii) (I — DDk) is invertible, and A is Hurwitz, where 



C K 
C ' 

Proof. Proposition 6 is a standard result [2, 21], though it is typically stated with the 
assumption that (C, A, B) is stabilizable and detectable. To prove the converse, note that 
if A is Hurwitz, then 

is stabilizable. 

It follows from the PBH test that (A, B) is stabilizable as well. A similar argument holds 
for detectability. ■ 

In this paper, we seek controllers that are both stabilizing and have a particular struc- 
ture. Therefore, we introduce new terminology to indicate this more restricted notion of 
stabilization. 

Definition 7. Suppose G G S(1Z P , k,m). We say that G is S -stabilizable if there exists 
K G S(1Z P , m, k) such that K stabilizes G. 

Note that these definitions are symmetric; K stabilizes G if and only if G stabilizes K. 
The same symmetry relationship holds for <S-stabilization. 

In the absence of structural constraints, it is well-known that every G G 7Z P can be 
stabilized. The main thrust of this paper is to explain what happens when G is structured, 
and we seek a controller that is 5-stabilizing. We will see in Sections 4-5 that not every 
G G S is iS-stabilizable, and 5-realizability plays an important role. 
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4 Simple Cases 

For certain systems, finding a structured realization is straightforward. In this section, 
we explore two such cases that will be useful later: stable systems, and linear chain 
structures. 
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4.1 Stable Systems 



If the plant is stable, it can always be iS-realized, regardless of the underlying graph. By 
duplicating states in a way that guarantees the correct structure, the resulting realization 
is always stabilizable and detectable since the starting plant was stable. We give the 
construction in the following lemma. A similar result appeared in [17]. 

Lemma 8. Suppose G G S(TZ'H <X> , k,m). Then G is S -realizable. 



Proof. We may construct a realization for G as follows. Separate G into its block-columns 
Gi for i = 1, . . . , N. Find minimal realizations 



G ? ; 





Bi ' 







In general, Ai and Bi will be full, but Ci and Di will have zero-rows corresponding to the 
zero- rows in d. Now stack the columns side-by-side and obtain a joint realization. 
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The realization (3) has the desired structure because A and B are block-diagonal, and C 
and D are in <S(M). Finally, the realization is stabilizable and detectable since we started 
with a stable system and each Ai is stable. ■ 

In the case where G is not stable, the construction method used in the proof of Lemma 8 
still produces a realization with the correct structure, but the realization may fail to be 
stabilizable and detectable. This is a consequence of the fact that realizing the block- 
columns separately and then re-combining them may cause some unstable poles to get 
duplicated. Indeed, this is precisely what happens if we attempt to realize (2). 



4.2 Chain Structures 



Linear chain structures correspond to adjacency matrices whose lower-triangular part is 
full. It turns out such systems are always 5-realizable. 

Lemma 9. Suppose G G S(Rp,k, m) and S has the full lower-triangular sparsity pattern 

\l 0" 



Then G is S-realizable. 

Proof. Let (A, B, C, D) be a minimal realization of G. We will sequentially construct a 
state transformation matrix T such that (T -1 AT, T~ X B, CT, D) has the desired sparsity 
pattern. Consider the simplest case, N — 2. Partition the realization according to the 
block-structure 
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G22 
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(4) 
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Note that D must already have the correct sparsity pattern. A realization for the zero- 
block G\2 is 



= 



A 


B 2 " 


Ci 






(5) 



Let T\ be the transformation matrix that puts (5) into Kalman canonical form. There 
are typically four blocks in such a decomposition, but there will only be three in this case 
since the system we are realizing is identically zero and thus there can be no modes that 
are both controllable and observable. Applying T\ to (4), we obtain 
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This realization is block-lower-triangular, and we notice that the index sets n may not 
be unique. For example, the modes A ss can be part of either the An block or the A 22 
block. 

For the case where N > 2, put the first block of the diagonal into Gu and lump the 
rest of the block- lower-triangular structure into G*22- Upon applying the T\ found from 
the N = 2 case, we are left with 



G 2 



Aqo 


Bcs 


C*23 


D22 . 



Apply this approach recursively by finding a T2 that puts the first zero-block-row of G22 
into Kalman canonical form. Continuing in this manner eventually yields a realization 
for G that has the desired sparsity pattern. Furthermore, this realization is stabilizable 
and detectable since it is minimal. ■ 



5 Main Results 

Our main results draw the connection between 5-stabilizability and 5-realizability. We 
will show that <S-stabilizable plants are always tS-realizable, and 5-stabilizing controllers 
are always 5-realizable. In other words, if structured systems connected in feedback 
results in a an internally stable interconnection, then both systems must be 5-realizable. 

The main result has two essential ingredients. In Section 5.1, we assume the plant is S- 
realizable. We then find a Youla-like parameterization of all 5-stabilizing controllers and 
show that all 5-stabilizing controllers are 5-realizable. In Section 5.2, we assume instead 
that our plant is 5-stabilizable, but make no assumptions regarding iS-realizability. We 
then show that there must exist a 5-stabilizing controller that is also <S-realizable. These 
results are combined into one concise statement, which we give in Corollary 14. 

5.1 Parameterization of Stabilizing Controllers 

In this subsection, we assume the plant is <S-realizable, and we give a characterization 
of all structured stabilizing controllers. This parameterization is similar to the classical 
work of Youla [19], and the unstructured version is one of the pillars of classical %2/Hoo 
theory [21]. In the framework adopted herein, the sparsity constraint imposed on the 
controller amounts to the same sparsity constraint being imposed on the Youla parame- 
ter Q. 
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Similar parameterizations have appeared in the literature. Most notably, the work of 
Qi et.al. [7] treat the lower-triangular case and other cases of interest. A more general 
parameterization is provided in [11], where the authors consider general quadratically 
invariant systems in which the plant and controller may have different sparsity patterns. 

Theorem 10. Suppose G G S(R p ,k,m) is S -realizable, with a structured realization 
given by (A, B, C, D). 

(i) G is S-stabilizable if and only if (Cn, An, Bn) is stabilizable and detectable for all 
i efl n . 

(ii) In this case, for i G Q n , choose F^ and Lj such that An + B^Fi and An + L^Cn 
are Hurwitz. Define F d = diag{Fi} ie o n and L d — diag{ij} ig ^„ • A particular 
S-stabilizing controller is given by 



A + BF d + L d C + L d DF d 


-L d ' 








Proof. Suppose (Cn, An, Bn) is stabilizable and detectable for all iGll. By definition, 
there must exist Fj and Li such that An + BnFi and An + LiCn are Hurwitz. Define 
F d = diag{Fi}i £ fj„ and L d = di&g{Li}i e n n . Since A is block-lower-triangular, A + BF d 
and A + L d C are also Hurwitz. It is straightforward to check that K defined in the 
theorem statement above is a stabilizing controller with the correct structure. 

Conversely, suppose G is stabilizable by some K G S(1Z P , m, k). Let S be the adjacency 
matrix with a full lower-triangular sparsity pattern. Then K G <S(1Z P , m, k). By Lemma 9, 
K is 5-realizable. Let (A K , B K , Ck, D k ) be a stabilizable and detectable realization for 
K structured according to S. By Proposition 6, (/ — DD^) must be invertible, and A 
must be Hurwitz. Rearrange the columns and rows of A by taking the first sub-blocks 
of each large block, then the second sub-blocks, and so on. The resulting matrix is 
block-lower-triangular, and has the same eigenvalues as A. The main diagonal blocks are 



An 
A Kli 



+ 



Bn 
Bkh 
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-Da 
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Cku 




Cii 






for 



Applying Proposition 6 once more, we conclude that (Cn, An, Bn) is stabilizable and 
detectable. ■ 

Theorem 10 establishes a necessary and sufficient condition for the existence of a 5- 
stabilizing controllers for for iS-realizable plants. Our next result is a parameterization of 
all 5-stabilizing controllers. 

Theorem 11. Suppose G G S(lZ p ,k,m) is S -realizable, with a structured realization 
given by (A,B,C,D). Let 





A + BF d + L d C + L d DF d 


-L d 


B + L d D ' 


J = 


F d 





I 




— (C + DF d ) 


/ 


-D 



where F d and L d are defined as in Theorem 10. Also, define 
Fi(J,Q) = Ju + J 12 Q{I - JiiQY x Jii. 

(i) Every K G S(JZ p ,m,k) that stabilizes G is parameterized by K = Fi(J,Q), where 
Q G SifRMoa, m, k) such that I + DQ(oo) is nonsingular. 



(ii) Every stabilizing controller is S-realizable. 
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Proof. By the proof of Theorem 10, A + BF d and A + L d C are Hurwitz. The classical 
parameterization of all stabilizing controllers [21] is given by Fi(J, Q) where J is defined 
in the theorem statement above, and Q E IZHoo such that I + DQ(oo) is nonsingular. It 
remains to show that K E S if and only if Q E S. This is a straightforward consequence 
of Remark 3, and the fact that all subblocks of J are in S. 

We will now show that every stabilizing controller is 5-realizable. Note that Q E 
S(TZHoo,m,k), so we may apply Lemma 8. Let (Aq, Bq,Cq, Dq) be a structured real- 
ization for Q that is stabilizable and detectable. The associated controller has a realization 





' 4 11 


il2 
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Fi(J,Q) = 


Mi 


i 2 2 


B 2 




c\ 
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b 



which can be computed using the Redheffer Star-Product [21, §10]. For example, the A 
and B terms are: 

in = A + BF d + L d C + L d DF d - (B + L d D)(I + D Q D)-\C + DF d ) 
ii 2 = (B + L d D)(I + DqD^Cq 
Mi = -B Q (I + DD Q )-\C + DF d ) 
i 22 = A Q - B Q {I + DD Q )- l DC Q 

Bi - -L d + (B + L d D)(I + DqD^Dq 

B 2 = B Q (I + DDq)- 1 

Using Remark 3 and the fact that the realization for Q is structured, it is straightforward 
to show that the iy and Bi are each structured according to S. The same is true of the 
C\ and D. Stabilizability can be verified explicitly. For example, one can check that 



"ill 


il 2 " 
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'Bi~ 


i21 


A 22 _ 


A 



[C + DF d DCq] 



A + BF d BC Q 

Aq 



which is clearly Hurwitz. A similar argument applies for detectability. Finally, we can 
permute the states as we did in the proof of Theorem 11 to obtain a realization for Fi(J, Q) 
for which all matrices are in S. This shows that the controller is 5-realizable, as required. 



5.2 Stabilizable Systems are Realizable 

In this subsection, we assume the plant is 5-stabilizable, but make no assumptions regard- 
ing 5-realizability. In this case, we show that one can always construct an 5-realizable 
controller. 

Theorem 12. Suppose G E S(1Z P , k,m) is stabilized by some K E S(1Z p ,m, k). Then, 

(i) G is stabilized by K d = diag(i^n, . . . , Kmn)- 

(ii) K d is S -realizable. 

Proof. The proof is similar to the proof of Theorem 10. G and K may not be <S- 
realizable, so we embed their sparsity patterns in S D S, the set of block-lower-triangular 
transfer matrices. Apply Lemma 9 to find triangular realizations, and the result follows 
from Proposition 6. ■ 
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Theorem 12 provides a simple test for S-stabilizability; stabilizing each subsystem 
individually must also globally stabilize the system. Otherwise, no ^-stabilizing controller 
exists. 



Corollary 13. Suppose G E S(lZ p ,k,m 
controller that stabilizes Gu. Define Kd 
if and only if Kd stabilizes G. 



For i = 1, 
diag(ifi, . . 



..,7V, 
,K N ). 



let K t 6 K™* xk * be any 
Then G is S-stabilizable 



Proof. If Kd stabilizes G, then clearly G is S-stabilizable, since Kd € S. The converse 
follows from Theorem 12. ■ 

The main result of this paper is that stabilizable plants are always S- realizable. We 
state it below as a Corollary. 

Corollary 14 (Main Result). Suppose G € S(1Z p ,k,m) is S-stabilizable. Then G is 
S -realizable, and every K G S{1Z P , m, k) that stabilizes G is S -realizable. 

Proof. By Theorem 12, every S-stabilizablc plant G can be stabilized by an 5-realizablc 
controller Kd- Interchanging the roles of the plant and controller, Kd is an iS-realizable 
plant stabilized by G. By Theorem 11, G must be 5-realizable. Applying Theorem 11 
once more to G, we conclude that all stabilizing controllers must be 5-realizable. ■ 

Note that the converse of Corollary 14 is generally false. One can construct an S- 
rcalizable plant that is not tS-stabilizablc. Consider for example: 
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If we realize G by clustering the unstable mode with An, i.e. using the index set n = (2, 1), 
it leads to an undetectable (Cn,An). If we instead associate the unstable mode with 
A 2 2 by using n = (1,2), it will lead to an unstabilizable (^22,^22). By Theorem 10, 
this plant is not stabilizable by a lower-triangular controller. Of course, non-structured 
stabilizing controllers exist since the realization for G given above is minimal. 



6 Concluding Remarks 

The main contribution of this paper is to establish the connection between the structured 
versions of stabilizability and realizability. In the context of structured controller syn- 
thesis, it is a basic requirement that the plant is 5-stabilizable and that the controller 
is iS-stabilizing. By Corollary 14, it is necessary that the plant and controller are both 
iS-realizable. The counterexample (2) causes no problems because as a plant, Gi would 
not be iS-stabilizable, and as a controller, G2 could not stabilize any structured plant. 

As mentioned in Section 2, the work of Vamsi and Elia [17, 16] also addresses the 
issue of structured realizability. The underlying assumptions arc different however. The 
authors assume discrete-time subsystems with one-timestep delays along directed edges. 
In the work herein, we assumed a delay-free network of continuous-time systems. Some 
results, such as the counterexample (2), are valid in either framework. However, the other 
results in this paper do not immediately translate to discrete-time. 

In this paper, we adopted the classical definition of internal stability. This notion of 
stabilization is the weakest one possible, since it does not guarantee that the signals that 
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travel between subsystems or between subcontrollers remains bounded. It follows that 
Corollary 14 will remain true even for stricter notions of stabilization. The work of Yadav 
et.al. [18] explores other types of stability in great detail. 

While 5-stabilizability guarantees 5-realizability, the question of how to explicitly con- 
struct a structured realization is still open. Lemmas 8 and 9 provide constructions for 
the special cases of stable and triangular systems respectively, but no universal method 
for constructing structured realizations is known. 
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A Proof of Counterexample 

In this section, we provide a proof that the example (2) of Figure 2 is not ^-realizable. 
Proof. Suppose G2 has a stabilizable and detectable realization in <S 2 - Label the blocks 
of this realization as 



A22 
A32 
A i2 



'A n 


A31 
An 






-433 








A44 



and similarly for B and C. Expand C(sl— A) 1 B and equate with G2. The four equations 
corresponding to each of the nonzero entries Gij are 



[Qj Cu] I si - 



33 







B 

Bj 



33 



(6) 



for (i, j) G {(3, 1), (3, 2), (4, 1), (4, 2)}. A minimal realization for G2 is given by 





" 1 


1 


1 

























G 2 = 



















1 
















1 















Since we have assumed our realization to be stabilizable and detectable, the eigenvalues 
of A must consist of the unstable eigenvalue 1, together with some number of stable eigen- 
values. Since A is block-lower-triangular, the eigenvalue 1 must appear as an eigenvalue 
of exactly one of the An for i = 1, ... ,4. For = (3, 1), we deduce from (6) that the 
eigenvalue 1 is contained in either An or A 33 , but for = (4, 2), we deduce that it is 
contained in either A22 or A44. This contradiction implies that no structured stabilizable 
and detectable realization of G2 can exist. ■ 
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